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We report on results of an experimental search for the Kl — ► 7r + 7r _ e + e _ decay 
D | mode. We found 13.5 ± 4.0 events and determined its branching ratio to be (4.4 ± 

1.3(stat.)±0.5(syst.)) x 10~ 7 . The result agrees well with the theoretical prediction. 
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1 Introduction 



In the previous Letter [1], we reported the results of our experimental search 
for the decay mode — > 7r + 7r _ e + e _ , which established the upper limit of 
4.6 x (90% CL) on its branching ratio. In this decay mode, CP violation 
may occur as an interference between two intermediate states with different CP 
properties [2,3] . Thus it can provide a new testing ground for investigating CP 
violation. A theoretical model predicts a branching ratio of about 3 x 1CT 7 [2]. 
Encouraged by this value, we continued the experiment further aiming to 
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establish this decay mode. Recently a group at the Fermilab (KTeV) reported 
the measurement of its branching ratio based on 36.6 ±6.8 observed events [4]. 
In this article, we present our new result [5]. 

The experiment was conducted with a 12-GeV proton synchrotron at High 
Energy Accelerator Research Organization (KEK). The experimental set-up, 
shown in Fig. 1, was already described in Ref. [1]. We used the same set-up 
with one vital change: the decay volume, which had been filled with helium 
gas in the preceding run, was evacuated down to 8 x 1CT 3 Torr. As will be 
seen, this was very effective to reduce backgrounds originating from nuclear 
interactions. Below we briefly describe our set-up for convenience. The Kl 




Fig. 1. Schematic plan view of the KEK-E162 detector 

beam was produced by focusing 12-GeV/ c primary protons onto a 60-mm-long 
copper target. The divergence of neutral beam was ±4 mrad horizontally and 
±20 mrad vertically. Following the 4-m-long decay volume was a charged parti- 
cle spectrometer. It consisted of four sets of drift chambers (DC1-DC4) and an 
analyzing magnet with an average horizontal momentum kick of 136 MeV/c. 
The drift chamber position resolution was in average 210 jum, and its effi- 
ciency was greater than 98% under a typical running condition. A threshold 
Cherenkov counter (GC) with pure N 2 gas at 1 atm was placed inside of the 
magnet gap to identify electrons. For the present decay mode, we obtained, by 
adjusting software cuts in the off-line analysis, an average electron efficiency 
of 99.7% with a pion- rejection factor of 50. Two banks of pure Csl electro- 
magnetic calorimeters, located at the far end, measured energy and position 
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of electrons and photons. Its energy and position resolutions were found to be 
approximately 3% and 7 mm for 1-GeV electrons, respectively. There were 4 
sets of trigger scintillation counters, called TC0X,TC1X,TC2X/2Y and TC3X, 
where X(Y) represented a horizontally(vertically)-segmented hodoscope. The 
trigger for the present mode was designed to select events with at least three 
charged tracks which included at least two electrons. 



2 General Analysis and Normalization Mode 

In order to exploit full statistical power of the available data, we combined the 
present data (vacuum data) with the one presented in Ref. [1] (He data), and 
analyzed them without distinction. Approximate Kl flux ratio of the vacuum 
to He data set was 5:1. We closely followed the off-line analysis procedure 
described in Ref. [1] with refinements at several places. At first, events were 
required to have 4 (and only 4) tracks with a common vertex in the beam 
region inside the decay volume. Then particle species were determined. A pion 
was identified as a track which could project onto a cluster in the calorimeter 
(a matched track) with E/p < 0.7, where E and p were an energy deposit 
measured by the calorimeter and a momentum of the track determined by 
the spectrometer, respectively. An electron was identified as a matched track 
with E > 200 MeV, 0.9 < E/p < 1.1, and GC hits in the corresponding 
cells. The events containing n + , e + and tracks were selected, and then 
applied two kinematical cuts. One was the requirement for the invariant mass 
of e + e~ pair (M ee ) to be at least 4 MeV/c 2 ; this cut was effective to suppress 
backgrounds from external conversion of 7-rays into e + e~ pair. The other was 
a limit on the pion momentum asymmetry (A + _ = (p 7T + — p 7r -)/(p n + + Pn-)) 
within ±0.5. This cut was introduced in Ref. [1] to suppress backgrounds due 
to nuclear interactions, but was found also useful to remove events with pion 
decays. The remaining events are our basic event samples for further analysis, 
and mainly contain both signal mode Kl — > 7r + 7r~e + e~ and normalization 
mode K L — > 7r + 7r~7r (n — > e + e~7). 

We now identify the normalization mode (denoted as 7r + 7r~7r£ in the follow- 
ing). We required at least one 7-ray with energy above 200 MeV. Here a 
7-ray was defined as a cluster in the calorimeter which did not match with 
any charged tracks. If more than one 7-ray existed, we selected the one for 
which the invariant mass of e + e~7 (M ee7 ) was closest to the tt° mass. Then 
we imposed three major kinematical cuts to select 7r + 7r~7r£, events. First, the 
invariant mass of e + e~7 (M ee7 ) was required to be within 3(7^ °f the n ° m ass 
(Mn.0), where the 7r° mass resolution om „ was measured to be 4.6 MeV/c 2 . 
Second, 9 2 was required to be less than 20 mrad 2 , where denotes the an- 
gle of the reconstructed Kl momentum with respect to the line connecting 
the production target and decay vertex. Fig. 2 shows the invariant mass dis- 
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Fig. 2. The invariant mass distribution of ir + ir e + e 7. The arrows indicate the 
region to select the Kl — > t[ + t:~tt® d events. 
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for the events which satisfied all the cuts 



mentioned so far. A clear Kl peak can be seen with a mass resolution (<Jm Kl ) 
of 5.5 MeV/c 2 . As a final cut, events were requested to lie within 2>(Tm Kl of the 
K L mass (M Kl ). After all the cuts, 12212 events remained [6]. We call these 
events the n + n~n° D reconstructed events. 



3 Signal Mode and Background Subtraction 



For the signal mode Kl — > 7r + 7r~e + e~, the major backgrounds came from 
7r + 7r - 7T£j events with an extra 7-ray missed from detection. We employed the 
same strategy as in Ref. [1]. At first, we removed events which had extra 
7-ray(s) consistent with n° — > e + e~7. To reduce the data set size, we also 
applied at this stage loose kinematical cuts; 9 2 <100 mrad 2 and 410 < M nnee < 
590 MeV/c 2 , where M n7ree being the invariant mass of 7r + 7T~e + e~. Then we 
defined a parameter which quantified consistency of an event with 7r + 7r~7r5). 
The parameter, called Xdj was gi ven by 




where p 7 was an arbitrary momentum of a 7-ray assumed to exist [7]. The 
standard deviations Om , &m k and a e [8] were all obtained from the 7r + 7r _ 7r^ 
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reconstructed events. We determined p 7 by minimizing xb- Fig. 3 shows the 




Fig. 3. The distribution of the square root of xb- The solid and dashed lines 
are the results of Monte Carlo (MC) simulation, respectively, for the tt^^iTq and 
7r + 7r~e + e~ modes. The Kl flux in the simulations was normalized with the recon- 
structed it + it~it° D events, and the branching ratio for 7r + 7r~e + e~ was assumed to 
be given by our final result. The arrow shows the cut position. 

\Jxb distribution of the events after minimization. (Here only the events with 
6 2 <20 mrad 2 and \M nnee — M^ L \ < 5cm Kl are plotted to illustrate clearly 
existence of two components.) The solid and dashed lines in the figure show 
the results of Monte Carlo (MC) simulations, respectively, for the it + it~it% 
and 7r + 7r~e + e~ modes. The Kl flux in the simulations was normalized with 
the reconstructed n + 7i~7r ( b events, and the branching ratio for 7r + 7r~e + e~ was 
assumed to be given by our final result. We considered the events with xb < 17 
as 7T + 7T~7r% events, and removed from the sample. The cut rejected 99.7% of 
the n + 7i~7i < b background, while retaining 92% of the signal. Fig. 4(a) shows 
the M^nee vs $ 2 scatter plot of the 7r + 7r~e + e - candidate events after the xb 
cut. Here the solid dots represent the vacuum data and the plus the He data. 
Let us first consider the vacuum data. We notice that there exists a cluster 
of events (15 events) inside our signal box defined by 9 2 <20 mrad 2 and 
~ M Kl \ < 3 cm Kl - Fig. 4(b) is the corresponding plot obtained by the 
ir+ir^ir'b MC simulation. Note that the MC statistics in the plot is about 5 
times more than the vacuum data. We also notice that there are non-negligible 
backgrounds remaining similarly in both plots. They are more or less uniformly 
distributed along the # 2 -direction, and form two distinct bands. It is found 
from the MC study that there exist two mechanisms for the n + 'K~7i < b events 
to pass through the xb cu ^- O ne i s 7-radiation: events in which e + /e~ radiates 
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Fig. 4. (a) The M n7ree vs 6 2 scatter plot of the Kl — > 7r + 7r~e + e~ candidate events 
after the Xd cut- The solid dots represent the vacuum data and the plus the He data. 
The box indicates the signal region, (b) The corresponding scatter plot obtained by 
the 7r + 7r - 7r£, MC simulation. Note that the MC statistics is 5 times more than the 
vacuum data. 



7-rays internally or externally would loose its energy. Such events constitute 
the low mass background band. The other is 7T-decay: events in which 
decays inside the magnet would have wrong momentum assignment. The high 
mass backgrounds are due mainly to these events. Next we consider the He 
data. There is no event left inside the signal box [9], as was reported in Ref. [1]. 
The backgrounds outside the signal box are spread more broadly compared 
with the vacuum data. These events are most probably due to the nuclear 
interactions off the helium nucleus. 



We now want to estimate the number of backgrounds inside the signal box. 
To this end, we relied upon the fact that, for both He and vacuum data, the 
background distribution was approximately flat along the # 2 -axis. We thus 
projected the events with \M W7Tee — M Kl \ < 3ctm Kl onto the # 2 -axis. The 
resultant distribution is shown in Fig. 5. There is a clear signal peak at 9 2 = 0. 
The solid line in the figure represents the 7r + 7r _ 7r£) MC simulation, in which the 
Kl flux was normalized with the reconstructed 7r + 7r~7r^) events. We denned 
a "control" region to be 30 < 9 2 < 100 mrad 2 . We then determined a scale 
factor by normalizing the number of 7i + n~i{ D MC events in this region (5.2 
events) to that of the data (6 events). The number of backgrounds was given 
by the number of -K + 7i^7i% MC events in the signal region (1.3 events) times 
the scale factor (1.15). We found this to be 1.5 ± 1.0, where the error includes 
uncertainty due to the MC statistics. Subtracting this from the events inside 
the signal box, we determined the number of signal events to be 13.5 ± 4.0. 
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Fig. 5. The 2 projection of the Kl — > 7r + 7r~e + e~ candidate events. The solid line 
in the figure represents the ir + ir~7r^ MC simulation. Its Kl flux was normalized 
by the reconstructed it + it~it® d events. The arrows indicate the signal and control 
region boundaries. 



4 Branching Ratio and Systematic Errors 



The branching ratio was calculated by 



Br(K L -> n + n-e + e-) = Br(K L -> Tr+Tr^Tr ) x Br(n° -> e + e _ 7) 
A(7r + 7r _ 7r£,) 77(7r + 7r _ 7r£,) iV(7r + 7r~e + e~) 
74(7r + 7r~e + e~) ?](7r + 7r _ e + e _ ) iV(7r+7r~7r£)) 

where A, 77 and iV denote acceptance, efficiency and observed number of events, 
respectively. The detector acceptances were determined by MC simulations, 
and were found to be 0.98 x 1CT 3 for 7r + 7r~7T£ and 2.6 x 1CT 3 for 7r + 7r~e + e~ [10]. 
Here, the latter acceptance depends upon the 7i + n~e + e~ matrix elements; 
we employed a theoretical model given by Ref. [2]. Most of efficiencies were 
common to the both modes, and they tended to cancel out in the efficiency 
ratio. The largest difference stemmed from the detection efficiency of extra 
7-rays in 7i + ii^7T^. Using a MC simulation, we estimated its inefficiency to be 
20.4%. The Xd cu ^ an d n ° ~^ e + e~7 inclusive cut were applied only to the 
signal mode: the efficiency of the former was found to be 92% while the over- 
veto due to the latter was estimated to be 1%. We found that all other effects 
caused negligibly small difference in the ratio of efficiency. Using the known 
branching ratio of Br{Ki — > 7r + 7r~7r°) = 0.1256 and Br(-K° — > e + e~7) = 
1.198 x 10~ 2 [11], we finally arrived at the branching ratio of 

Br(K L -> 7i + n-e + e-) = (4.4 ± 1.3 ± 0.5) x 10~ 7 , 
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Table 1 

Summary of systematic errors in the branching ratio. 



Source 


uncertainty 


Kl momentum spectrum 


4.8% 


Matrix element 


3.9% 


Others 


3.1% 


Background subtraction 


7.4% 


Nuclear interaction 


3.6% 


Other contamination 


1.4% 


Br(K L -» vr+Tr-vr^) 


3.1% 


Total 


11.3% 



where the first (second) error represents the statistical (systematic) uncer- 
tainty. 

Table 1 summarizes our systematic errors, which are divided into three cate- 
gories. The first one is related to the acceptance-efficiency ratio. The largest 
contribution in this category came from uncertainty in the Kl momentum 
spectrum employed as an input to the various MC studies. It affected both 
acceptance and efficiency, and resulted in the fractional uncertainty of 4.8% 
in the final branching ratio. The second largest contribution stemmed from 
uncertainty in the matrix element of the theoretical model. This brought 3.9% 
uncertainty in the acceptance ratio. The second category is related to the num- 
ber of the events iV(7r + 7r~e + e~) and/or N(tv + tv~tv l \). We checked stability of 
the background subtraction by employing different control regions. Actually, 
widening the region in the M n7ree direction, we tried \M n7ree — Mk l \ < ^><Jm Kl 
and < 7<7m Kl instead of < 3cxm Kl - We considered the biggest difference re- 
sulted from this treatment as a systematic error (7.4%). For the He data, we 
also tested different event shapes assumed for the backgrounds due to the nu- 
clear interaction. The effect was at most 3.6% as long as they gave statistically 
consistent results with the null observation. All other background contamina- 
tions, such as external 7-ray conversion events associated with Kl — > 7r + 7r~7 
and K L — > 7r + 7r~7r° (n° — > 77), were estimated to be small (1.4%). The final 
category is related to the current experimental errors on the branching ratios 
Br{K L — > 7T + 7r~7r°) and Br(n° — > e + e~^). Summing up all the uncertainties 
in quadrature, we found the over all systematic error to be 11.3%. 

In summary, we observed 13.5 ±4.0 Kl — > 7r + 7r~e + e _ events, and determined 
its branching ratio to be (4.4±1.3(stat.)±0.5(syst.)) x 1CT 7 . The result is found 
to agree with the theoretical prediction as well as the recent measurement. 
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